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Abstract 

The electromagnetic theory of the strongly driven ion-temperature-gradient (ITG) instability 
in magnetically confined toroidal plasmas is developed. Stabilizing and destabilizing effects are 
identified, and a critical j3 e (the ratio of the electron to magnetic pressure) for stabilization of the 
toroidal branch of the mode is calculated for magnetic equilibria independent of the coordinate along 
the magnetic field. Its scaling is /3 e ~ Lxe/R > where Lx e is the characteristic electron temperature 
gradient length, and R the major radius of the torus. We conjecture that a fast particle population 
can cause a similar stabilization due to its contribution to the equilibrium pressure gradient. For 
sheared equilibria, the boundary of marginal stability of the electromagnetic correction to the 
electrostatic mode is also given. For a general magnetic equilibrium, we find a critical length (for 
electromagnetic stabilization) of the extent of the unfavourable curvature along the magnetic field. 
This is a decreasing function of the local magnetic shear. 
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I. INTRODUCTION 


Most kinetic investigations of ion-temperature-gradient (ITG) instabilities in plasmas 


notable exceptions [f 
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a numerical standpoint 


rely on the simplifying assumption that perturbations are electrostatic BA. With some 
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ectromagnetic perturbations have been considered mainly from 
I 5 J], and the attempt to understand their role in ITG stability 


has resulted in a patchy collection of numerical findings rather than in a coherent physical 
picture. Moreover, most studies have neglected magnetic compressibility and have thus 
neglected the magnetic perturbations parallel to the equilibrium magnetic field, SBn, that 
are generated by the instability to maintain perpendicular pressure balance. 

From an analytical point of view, the equations describing electromagnetic ITG modes 


are present in the works of Antonsen and Lane 4], and Tang, Connor and Hastie jfjj], where 
the linear theory of kinetic ballooning modes was formulated. However, since the equations 
derived in these papers are general and thus encompass many types of instabilities, the role 
of ITG modes is somewhat obscured. Later, Kim et al. Jfi| focused on the physics of the 
toroidal ITG instability, extending previous electrostatic work to finite f3 (the ratio of the ion 
to magnetic pressure) by including the effect of the induced electric held, —dA\\/dt, on the 
electron motion along the magnetic held. Therefore, effects of A\\ (the parallel component 
of the magnetic vector potential) were included, however, effects of 5B\\ were neglected on 
the grounds that f3 was considered to be small. 

A complete electromagnetic theory of ITG modes must retain all three gyrokinetic helds: 
0 (the electrostatic potential), An and SB y. Formally, the latter two are finite-/3 effects, 
but as we shall see they become important at surprisingly low values of (3 because of other 
small parameters present in the problem. A general finite-/3 theory must necessarily describe 
several families of instabilities, such as ITG Alfvenic modes |7|, 0—induced Alfvenic eigen- 


modes [16], /?—induced temperature gradient eigenmodes [13], and kinetic [4 5] and ideal 
ballooning modes jl8|. In this work, we limit ourselves to the analysis of curvature-driven 
ITG modes by adopting an ordering scheme which excludes other instabilities but, at the 
same time, allows a small value of [3 {[3 -C 1) to affect the ITG mode through both An and 
SB\\. The result is a simple formulation shedding light on why and when electromagnetic 
effects are important for toroidal ITG instabilities. 

From a numerical point of view, early gyrokinetic simulations ll|, [12] had already found 
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magnetic compressibility to be important, in particular to cancel the stabilizing effect of 


the “self-dug” magnetic well for drift instabilities 


19]. Waltz and Miller reported on such a 


cancellation, resulting in a substitution rule for the magnetic drift: magnetic compressibility 
could be dropped if the magnetic drift were replaced by the curvature drift jll|. While this 
fact now seems to be common knowledge in part of the gyrokinetic community 14, |20j, the 
picture that emerges from systematic electromagnetic gyrokinetic simulations of the ITG 


13] and difficult to disentangle. A simple analytical explanation 


mode is more complicated 
is therefore helpful. 

In the present work, we build on the recent electrostatic linear theory of Plunk et al. 


2l|, exploit asymptotic techniques to solve the kinetic problem of the ITG instability, and 


identify the conditions that allow this theory to accommodate electromagnetic perturbations. 
Somewhat to our surprise, we find that, for strongly driven modes, magnetic compressibility 
can be as important as perpendicular magnetic perturbations for values of /3 accessible to 
both tokamaks and stellarators. The ions contribute to magnetic perturbations to maintain 
pressure balance, whereas the electrons can have both a stabilizing and destabilizing effect, 
depending on the value of /3. In the case of a uniform equilibrium magnetic field, a new 
critical (5 for the electromagnetic stabilization of the toroidal ITG is calculated. This differs 
from the one given by Kim et al. 6] in a fundamental way. A similar stabilization is 
predicted when an additional fast particle population is considered. For sheared magnetic 
equilibria, the boundary of marginal stability for the electromagnetic component of the ITG 
is given, for the first time, using a local approximation of the magnetic drifts. 


II. PHYSICAL PICTURE 


To understand the role of magnetic perturbations for ITG modes, it is useful to start with 
a physical picture of the instability. We follow Rosenbluth and Longmire, who first described 


the physical mechanism responsible for interchange modes 


22]. The same description works 


for the curvature-driven branch of the ITG modes and will be used here. 

Consider a plasma with gradients of the temperature and the magnetic field strength in 
the direction of —Vx. The magnetic field points in the ^-direction, and for simplicity we 
take the density gradient to vanish. The ion guiding centers drift in the direction B x VB, 
i.e., in the negative y-direction, and do so with a speed that decreases with increasing x, 
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since the drift velocity is proportional to the energy. 

If the plasma is displaced by an E x B drift in the x-direction by the distance 

£ = x£ 0 sin (k±y ), (1) 

the ion pressure is perturbed according to 

Spi = -£ ■ Vpi, 


where p t is the equilibrium ion plasma pressure. The ion guiding centers will then start 

accumulating at k±y = 2rm and a corresponding deficit of ion guiding centers arises at 

k\_y = (2n + l)7r, see Fig. (pQ). An electrostatic potential, 0 = 0 O cos k±y, thus builds up 

(with 0o having the same sign as £ 0 ) and gives rise to an E x B drift, 

d£ b x V0 0o. _ , , fc±0o . , k ± (f )o 

— = ——— = —z x Vcos k±y = x—— sm k^y = — 


that amplihes the initial perturbation (CP). In this picture of the instability, the motion of 
the ions parallel to the magnetic held is neglected, so it is tacitly assumed that k\\v t hi "C ca, 
where v t hi = ( 2Ti/rrii ) 1//2 denotes the ion thermal speed and co/k\\ the parallel phase velocity 
of the instability. The electrons, on the other hand, can be expected to move quickly 
compared with the instability, knv t h e w, and will therefore only experience a small E x B 
displacement. 

How is this mechanism affected by electromagnetic terms within the gyrokinetic descrip¬ 
tion of the instability? As already mentioned, there are two such terms, proportional to 
A\\ and 5B\\ (the perturbation of the magnetic held strength), respectively. The hrst one 
describes the effect of the inductive electric held and is important to the electrons, which 
unlike the ions have time to move significant distances along the magnetic held during the 
evolution of the instability. They are therefore sensitive to the parallel electric held, 

dAn 

B 'i = - V ii *--aT 

Instead of An, we introduce the quantity 0, dehned by 

8An 

ViV- = --gp 

so that An = —V|| (0 — 0). Ampere’s law, kj_An = poJ\\, then implies 

s v ' 7 » 
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where Vn J\\ describes the local accumulation of electrons due to their parallel motion. If the 
inductive field is weak, -0 <E </>, the electrons are approximately Boltzmann-distributed, 


8n 

n 


e( t>o , 

—- cos k±y, 
T ±y, 


so that 


Hence and from Eq. 


dSn ne 2 <9</> 0 

Vii Jii ~ e-^7- = ——— cos k±y. 


dt 


T P dt 


we obtain the estimate 

/done 2 d 2 (fto 


$ 


k 2 k 2 T P dt 2 


cos k±y. 


and we conclude that the critical (3 = 2 /iquT/B 2 above which electromagnetic effects are 
important, -0 ~ 0, scales as 


A 


k±Pik\\c s 


CO 


(3) 


where c s denotes the sound speed, and pi the ion Larmor radius. The ITG mode has 
a frequency of order cu* ~ (k±pi)c s /Lj_, where L± is the length scale of the cross-held 
gradients, it thus follows that the critical beta is j3 c ~ (k\\L ± ) 2 ~ e 2 , and, in a standard 
tokamak, can be ordered as the square of the inverse aspect ratio. This is the basic reason 
why electromagnetic effects are already important in standard tokamak situations when 
j3 ~ 10~ 2 rather than when j3 = 0(1). This critical beta also defines the value at which 
kinetic Alfven waves are relevant, since 


/3 C 


/ k±p s k\\v A 

V w 


/3, 


(4) 


and oo = k±p s k\\i>A is the kinetic Alfven wave dispersion relation. 

The other electromagnetic term in the gyrokinetic equation involves 5B\\ and is sometimes 
neglected in analytical treatments and numerical simulations of the gyrokinetic equation. 
Physically, it accounts for the perturbation in the VB drift due to the variation in magnetic 
held strength, 

V£> = V |B + <SB| ~ V{B + 6B\\). 


The latter is determined by perpendicular pressure balance, 

B 2 


<3 \p± + 


2p 0 


= 0, 
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with p± = p± : i + p_ l i6 , which implies 


SB I, = 


ho Sp± ho p'So 


B B 

and thus gives rise to a perturbed VB-drift of the ions 


sin k±y, 


= 


20,B 


b x V< 5 B|| ~ —x 


kovlp'i^o 
2 QiB 2 

’ v" 

/ 


cos 


(5) 


where b = B /B and O, = eB/rrii. As is clear from Fig. (JTJ) , this extra drift reinforces the 
density accumulation around k±y = 2mr and thus amplifies the instability. It does so even 
in the absence of a density gradient, since the VB-drift is proportional to the perpendicular 
kinetic energy and we assume that a temperature gradient is present. 

There is, however, also a third effect of finite plasma pressure, since this affects the 
equilibrium magnetic held by making the curvature vector deviate from the gradient of the 
held strength, 

_ Vj_B y 0 X7p 
K ~ ~ + 

where p = Pi + p e - The equilibrium ion drift velocity can thus be written 
,2 Hxk 


Vd = ( Y + V\ 


a 


ho^l 

2 0,B 2 


b x V ( pi + p e ) = v K + Av d 


( 6 ) 


— 

/ / 


hi 


IV 


where the second term on the right (III + IV) opposes the basic curvature drift (term 
II) causing the instability. Thus, if the plasma pressure is increased whilst the magnetic 
curvature is kept hxed, then the drift velocity is reduced and the instability is weakened. As 


has been discussed in the literature 


B, Q, 


14, |23|, this effect from the ions partly cancels 
that from 5By, but it is important to keep in mind that this cancellation only holds if n, 
rather than Vj_B, is held constant. A simple mathematical argument for the cancellation is 
given in an Appendix. 

Even though the electrons contribute relatively little to the ion instability, their pressure 
gradient exerts a stabilizing effect. While term III in Eq. ((HD tends to cancel the perturbed 
grad-B drift / of Eq. (J5j) , the diamagnetic electron contribution [term IV in Eq. (JH|) |, tends 
to oppose the drive of the mode //. When these terms balance, we have 

L r 


Pe 


J p 

R ’ 


(7) 
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where we used VB ~ B/R, and Vp e ~ po/L p . This stabilizing influence of finite [3 was 
studied by Hastie and Taylor for MHD instabilities in a combined mirror-cusp magnetic 
configuration [24], and by Rosenbluth and Sloan for electrostatic and weakly electromagnetic 
instabilities 25j. It will be confirmed quantitatively in the context of the electromagnetic 
ITG instability below. 

It is worth noticing that a similar stabilization can be expected when a population of fast 


ions is present 


26] . Just like the electrons, fast ions move quickly along the magnetic field, 


and in addition they have large gyroradii. If the typical velocity of the fast ions exceeds the 
phase velocity of the instability along the field, uj/k\\Vf ast < 1 , or their gyroradius exceeds 
the perpendicular wavelength, k±pf ast > 1 , such ions will experience relatively small ExB 
displacement and therefore contribute little to the magnetic-drift perturbation in Eq. ([5]). 
The fast ions will then contribute relatively little to the instability. On the other hand, their 
equilibrium pressure can be significant and acts to reduce the equilibrium drift in Eq. (jfJJ) by 
a new additive term giving pi + p e —* pi +p e + pf as t ■ We thus expect a net stabilising action 
from fast ions. Gyrokinetic simulations of plasmas with such particles indeed indicate the 
presence of a critical f3 for electromagnetic ITG stabilization that decreases with L p / R .27]. 
Moreover, the effect of fast ions is even more significant in nonlinear simulations. 


III. REDUCTION OF THE GYROKINETIC EQUATIONS 


Bearing in mind the qualitative picture from the preceding Section, we now give quan¬ 
titative substance to our findings. We proceed by first deriving from gyrokinetics a set of 
second order differential equations for the electrostatic and the magnetic potentials. These 
equations support the electrostatic ITG mode in the limit of vanishing /3, Alfvenic pertur¬ 
bations, magnetic compessibility and finite-ion-Larmor radius effects. They are derived in 
a large-7^ expansion, where rji = dlogR/dlogrii = L n . / Lr, , with T* and n* the equilibrium 
temperature and density, respectively. Kinetic ballooning modes are therefore diamagneti- 
cally stabilized within our ordering. 


Our starting point is the linearized gyrokinetic equation in ballooning space 


jd, 5, 


281 
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( 8 ) 


eF ( 


iv\\V\\h s + (u - Lu ds ) h s = (u - ul) — 

J-c 


Os 


T ( \ (x A \ , Ts n V ± J 1 (°*) 5B \ 

Jo ( a s) \<p ~ ^ll^ll) H-2—- 


J ths 


B 


where </> is the electrostatic potential, Ah the perturbed magnetic potential parallel to the 
equilibrium magnetic held in the Coulomb gauge, V • A = 0, 5B\\ the parallel magnetic 
held perturbation, and B the modulus of the equilibrium magnetic held. The form of the 
perturbations used is ~ exp [—icot + zkj_ • x]. The function h s , dehned by /r s exp (— iL s ) = 
Sf s + Z s e(j)F 0s /T s , denotes the nonadiabatic part of the perturbed distribution function, Sf s , 
where f s = Fq s + 5f s , with 5f s Fq s , Fq s is a Maxwellian equilibrium with temperature 
T s = rn s v't hs /2 and density n 0s , L s = k x vj_ • b/f2 s , b = B /B, with vj_ the perpendicular 
particle velocity. Here = Z s eB/m s is the cyclotron frequency, Jq and J\ are Bessel 
function of the hrst kind of argument a s = v±k±p s = v±y/2b, where p s = v t h s /kl s is the 
Larmor radius, k\ = k 2 (l + s 2 z 2 ), with k y the mode wave number, s the local magnetic shear 
and z the distance along the equilibrium held lines. Furthermore, Cj ds = 2 [u>bv\/2 + oj k v jj^ , 
2ub = kj_p s • Vthsb x V-B/B, 2uj k = kj_p s • v t hs b x ^b ■ Vbj , with uy the parallel particle 
velocity. Finally, uuj s = oj* s + p s uj* s (■ v 2 - 3/2), v = v/v ths , and uj* s = (1/2 )kyp s v ths /L ns . 

The gyrokinetic equation (JHJ) is most easily solved for the electrons, which we take to be 
sufficiently light that the terms multiplied by un dominate. Neglecting magnetic trapping, 
we thus obtain the electron response being described by the solution 


h P « - ( 1 * ^ 

UJ 


eijj 

T; 


Oej 


(9) 


where we have written Vii^ = iuA\\. 

For the ions, Eq. (JS]) is solved iteratively using the ordering 

k l v thi oj oj k + ojb , 

- - r\j - rs_/ Q r\j £ •C 

UJ 2 PiUJ*i UJ 



1, 


( 10 ) 


which retains the strongly driven (r/,; 1) toroidal and slab ITG instability and finite 

Larmor radius (FLR) effects. To include electromagnetic perturbations in the electrostatic 
picture, we use a maximal ordering for the helds, v t hiA\\ ~ e</>, and fold in lowest order 

h (o) = u-ujT 
UJ UJdi 
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Jo (a,) ^ + 2v‘ 

J-Oi 


Ji (a*) SB\\ 


Cli 


B 


on 


( 11 ) 













The electrostatic potential is obtained from the quasineutrality condition, 

n 0 e(T e _1 + T” 1 )^ = f d 3 vJ 0 hi - [ d 3 vh e , 


( 12 ) 


to which the contribution from /i;° becomes 


d 3 t-,/ n h (0) = 


e</> 


1 - (1 - riib) - 


u 




ViU*i $B\\ 
L0 B 


in lowest order. This density perturbation is a factor r /^ 1 -C 1 smaller than expected from 
the size of hf' ~ (^w*i/^)(e0/Tj)F O j, compelling us to find the solution to higher order. We 
thus iterate the solution, 


A 


U|| 


U|| 


(U - (U* 


LU — U7 


eAu 

i) do (a*) -^—Foi + zVii 


T 

- 1 l 


h 


(o) 


P & 

U - wJi) Jo (a*) ^7- 


+ * V ii h) 


( 0 ) 


(13) 


cu - cu* 

and find that a sufficiently accurate expression for the ion density perturbation is 


d 3 vJ 0 h, l = y 


i u\ ViU*i(uB + UJ K ) 

1-(1 - rub) - 


u 


U) z 


ViV*i SB\\ 

uj B rriiUJ 3 ^ 


eV||(0 — -0) 


(14) 


The magnetic field strength fluctuations are determined by the perpendicular Ampere’s 
law, 

= ^2 f d3 vm s v 2 ± a~ 1 J 1 (a 8 ) h s , (15) 

S 

implying that 5Bn/B is proportional to /3e0/Tj. At this point, a traditional and popular 
approach would be to neglect the magnetic compressibility altogether Q, 8-10, 15, 29|, since 
/3 in fusion relevant plasmas is of the order of 1% — 5%. However, even such a small j3 is 
not necessarily negligible, since it gets multiplied by a large factor of order e _1 in Eq. (T141) . 
In fact, using Eqs. (jUJ) and to calculate the integrals in Ampere’s law, we arrive at the 
conclusion that /3 ~ co 2 / (rjfujd) ~ e 2 1 is the correct ordering that allows us to calculate 


tJb 


(1) 


A ViU* 


1 Ve , 

h— ~4> 
T Vi 


e B 2 ' • ' r r ' ’ ( ' 16 ' ) 

where A = 2/i 0 nTi/B 2 and r = TjT e . This result is a special case of a general formula 
derived in the work of Tang et al. [5j on kinetic ballooning modes. Finally, using Eqs. (jUJ), 
(fl4l) and dT6l) in the quasineutrality condition (fT2lh we obtain 


r + 


_ A Ve rjjuti _ ViU*iVm d 2 
u 2 t r)i u> 2 2 o; 3 /m dz 2 


(0 - VO = 


2 rji - 5 - rji—b 

UJ Z L0 


0 (17) 
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where we have defined the normalization length L and the coordinate z along the held so 
that Z||V|| = d z . We have also used the result 

UJ K - UJ B = [1 + Ve/{Trji)]ViU*iPi/2, (18) 

which follows directly from the force balance equation Q 

j x B = Vp. (19) 


Equation (fI71) is similar to previous results in the literature, but is different in a couple of 
ways. In particular, the third term on the LHS is absent from previous electromagnetic 
theories of ITG instabilities j(]]. Another novelty of this equation is that the inclusion of 
the ion contribution to magnetic compressibility [the term proportional to 0 in Eq. 01611 ] 
resulted in the “rule” that the drive of the toroidal branch of the ITG [the first term on 


the RHS of Eq. (H 
numerical works 111 


is the curvature drift only. This result has been confirmed by various 

14, bo], 

to obtain 


We close the system of equations calculating the divergence of the current 


1 v m/ l \\ d ( d \ 

m~Tz 


00* 


1 7U 
T Pi 


( 20 ) 


This is obtained by using Ampere’s law after taking the T, s e s f d 3 vJo moment of the gyroki- 
netic equation (jBJ) 


B _ / k 2 , _ A ne 2 co^, , / 


IMP* 


> + £, 


OJ 


B 


d 3 vJ<y^-h S) 

OJ 


( 21 ) 


where the ordering in Eq. ([TO as been used. Velocity-space integrals are performed using 
solutions (0 and (ITT]) . Thus, from Eq. (12TT) . we have 


B 

PqOJ 2 


V 



PiU*i{uj K + 0J B ) 


OJ* 


CO*i i i 

Pi—bcp 
OJ 


„ 1 Pe 

1 H-- 


T Pi 

1 p e PiUJ*i(uj K + 0J B ) 


Vi¬ 


co * 


T SB, 

-L 7. 


(1) 


co e 


B 


r Pi 


OJ- 


A 


Substitution of Eqs. (fTO) and (fT 8 l) now yields Eq. (1201) . 


( 22 ) 
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IV. CRITICAL (3 FOR STABILITY 


Let us first consider the simple case in which the equilibrium magnetic held is independent 


of the coordinate along B. Then s = 0, kr^ = k 2 , and we can Fourier transform Eqs. (1171) - 
d2DD, to obtain 


r + 


cn* 


Pi r] e u>t , ojt k 2 z v 2 


UJ 


+ 


2 : u thi 


2t rji uj 2 2u uj 2 


A = 




uj z 


UJ 


where ojt = J]iUJ*i, 


and 


A = 


( 1 + (Pmhd - Pi) + 2 Z: b Pi 


JJe 

rrji 


. ( Pmhd Pi) T Pt, 


Pmhd — b- 


k 2 v 2 

V thi 


(23) 


(24) 


(25) 


2u K r]iUJ^[l + r] e /pn-]i)\ 

is the value of P above which ideal MHD modes would be destabilized if they were not 
suppressed by diamagnetic effects. 

For p MHD < Pi, A « rr]i[l + rie/prrji) - bu>/(2uj K )\/r] e , while for P MHD » ft, A « 1. 
Similarly, for Pi —» 0, A —> 1, and Eq. (f23l) reduces to the dispersion relation for the 

[ 13 , 0 , 


electrostatic ITG mode 


k 2 v 2 hi = ^ UJ T UJ K + b UJT 


T ^ -^ O 2 

u 2 uj uu z 


UJ z 


UJ 


(26) 


Equation (l23l) agrees with the large-r^ limit of Eq. (25) in Ref. 6] only if the electron 
contribution to the magnetic compressibility (the third term on the LHS) is neglected. In 
general, the coupling of all the roots of Eq. ([23]) is essential to understand the electromagnetic 
stabilization of the toroidal ITG mode. To illustrate a somewhat typical case, we solve Eq. 
(1231) numerically for r = 1, b = 0.1, R/L ^ = 5, k z l\\ = k z qR = 0.5, and q = \/2, where R is 
the major radius of the toroidal device, and q measures the pitch of the magnetic held. We 
consider the flat density limit for simplicity, uj ** = 0, but ujt p 1 0. For these values Pmhd = 
0.0125. We note the normalised frequencies ujt / (vm/qR) = y/b/2qR/L Ti = VbR/L T . and 
uj k / ( Vthi/qR) = Qy/b/2 = \Pb, for this particular value of q. Several electromagnetic branches 
can be observed, depending on the value of Q t . For small Pi we find two complex conjugated 
ion roots. In Fig. (]2]) we identify the toroidal ITG branch, 0 < 3ft[c o/{v t hi/qR)\ < 1, and 
0 < ^s[uj/( vthi/qR)\ < 1- Its P ~stabilization occurs at a critical Pf rit for which the imaginary 
parts of the two complex conjugated roots coalesce. At low /3*, a further stable electron 
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mode ^{uj/{vthi/qR)\ < 0 is present. Its real part changes sign when the f3 —stabilization of 
the ITG becomes effective for (3i ~ 1%, see Fig. (J2]). 

To establish the scaling of the observed P for stabilization with ojt, we solve Eq. (1231) 
for several values of lut / {v t hi/qR) = \]b/2qR/L T . but fixing u K /(v t hi/qR ) = QyJb/2 = \Pl, 
for the above values of b and q. This means ut / (vthi/qR) = V0.1 R/Lt,.. We then record 
the value of Pi at which the mode is completely stable. To determine the scaling of the 
observed P for stabilization with k z , we repeat the same evaluation of f3/ nt for constant 
ujt/( vthi/qR) = \Jb/2qR/L t { = V0H5, but varying k z Rq. As Figs. ([3])-(j4j) show, the critical 


(3 for stabilization scales as Pf nt P A 


PI 


, which implies 

1 l t . 


301 


crit 


(27) 


2q 2 R 

However, as is evident from the figures, Pf rtt lies somewhat above Pmhd-, which means that, 
for these parameters, the stabilization occurs only for values of Pi above the ideal MHD 
threshold. 

It is interesting to analyze the stability below this threshold, for /3i -C Pmhd- This situation 
corresponds to A « 1. For a strongly toroidal mode 4,oj t lo k k 2 v/ hi ujT / u, or 


Pi < Pt. 


<C b — « 
lot 





ujt 


(28) 


the new term on the LHS of Eq. (j23l) cannot be neglected, and indeed it is responsible for 
a new critical electron (3 e for stabilization. After neglecting the stabilizing FLR term on 
the RHS of Eq. ([23]) . we obtain tco 2 = —2lo t U) k /A + PiVe^r / (^ T Vi)■ Hence, the electron 
contribution to magnetic compressibility suppresses the instability when 

Vi 4 


Pe > PT U = 


-, for A > 0. 


(29) 


Ve A u T 

In the limit Pi Pmhd, A « 1. The same critical (3^ lt for stabilization is obtained in the 

Pmhd -C Pi limit, but now A ^ 1. In both cases, we End 


per it ^-‘T e 

A ~ ~R' 


(30) 


It is perhaps interesting to notice that f3 c i r " lt and P/ nt show different explicit scalings with 
oj k , however they follow the same scaling with R/Lj-- 

To verify the estimate in Eq. (1291) . we now solve Eq. (1231) numerically in the asymptotic 
regime ca*; = 0, k z qR = 0.001, ujt/( v t hi/qR) = 10, u> K /(v t hi/qR) = 0.25, r = 1, b = 0.05, and 
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q = 1.58. For these values A = 2—co /(2 v t hi/qR) ~ 2, when co /( vm/qR ) -C 1. Again, we solve 
Eq. (1251) for several values of ut/( vm/qR) at fixed co K /( v t hi/qR ) = 0.25, and co K /(ikhi/qR) = 
q^/b/2 at hxed cot/ (vm/qR) = 10, and record the value of fa at which the mode is completely 
stable. As Figs. (15|1 - (J6|1 show, the critical fa for stabilization agrees with Eq. (j29l) . In Fig. 
a, we notice a window of stability for (3 ~ 1%. A further destabilization might occur for 
/? > 2%. Incidentally, the high-/? unstable mode is not the kinetic balloning mode, since this 
is diamagnetically stabilized in our large-7/.,; limit. The presence of electromagnetic roots can 
be investigated further by considering the limit fa —>■ (3mhd■ 


A. Electromagnetic roots at fa « (3 MHD 


Equation (1231) . when co = 0, is in general a quartic for co. However, near fa ~ (3 MHD , it 
can be factored into a stable solution 


and a cubic 


with 


and 


^stab 


k Z V thi 


1 + 


Ve \ (3mhd fa 


1 + bW 


UJ CLiUJ + CLq — 0 , 


a 1 = — 


^6 + 2 ) 0 ;^, 


r yi + Tr/i/rje 4 u 2 K 

1 f uj T k 2 z v 2 hi , , u 2 u t 


ao — ~ 

T 


+ 4- 


2 b 

For fa —>• (3 MHD) co s tab —» 0, whereas the roots of Eq. (1321) are 


(31) 


(32) 

(33) 

(34) 


uj\ — A B , 

uj 2 = —- (A + B) + jl— (A- B ), 
and _ 

U3 = A(A + B)-i^-(A~B), 
with A 3 = C + VC 2 + D 3 , B 3 = C — VC 2 + D 3 , 

C — —-ao < o, 


(35) 

(36) 

(37) 

(38) 
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and 


For 


D = -ai < 0. 


f)V(f) 2 <o, 


(39) 


(40) 


all three roots are real. 

In the limit (ai/3 ) 3 <C (a 0 /2) 2 , A 3 ~ a 0 (af/ag) <C 1, and B 3 fa — a 0 = 0( 1). Thus, we 
find the unstable mode 


uj\ fa e*3 


1 ( uTkjvthi 
T 1 2 b 


1/3 


(41) 


while CU 2 = u\ is damped, and 013 = — |o;i| is marginally stable. In the case of negligible slab 
drive, we have 


oil fa e 




1/3 v m , Rb 2 ^27 

T J RWl'J 3 ° r L Ti r 2 ' 


(42) 


In the opposite limit (ai/3 ) 3 3> (ao/2) 2 , we have A 3 fa — i(|oi| /3) 3 / 2 , and T > 3 fa —A 3 . 
Therefore, we obtain one stable ion root 




1 k z V thi , 


_r \1 + rr]i/r) e 4a ; 2 
which, for negligible slab drive is 

kyPi Vth/i 


1 1/2 


6 + 2 co^coj 1 


(43) 




~ 1 1 -W k O > T 
T 


r R b 2 27 

2 ° r Z ^7 > T' 


The second root is 


CU 2 fa 0. 


(44) 


(45) 


Finally, we find the stable electron mode 


u 3 = -U>i- 


(46) 


V. ELECTROMAGNETIC BOUNDARY OF MARGINAL STABILITY. 

For the more realistic case of finite shear, we have k‘]_ = k 2 (1 + s 2 z 2 ) , and the previous 
analysis does not apply. Nevertheless, we can still construct a perturbative electromagnetic 
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if we use a local approx- 


theory of the ITG instability similar to that introduced in Ref. [^|, 
imation of the curvature drift, uj k (z) = u K (l — az 2 ) [21]. We calculate the electromagnetic 
correction to the electrostatic eigenvalue using a low [diUJ^/u 2 -C 1 subsidiary expansion. The 
zeroth order electrostatic response is given by Eq. (fT7D with = 0 , and q= exp[—Az 2 ], 
with [ 21 ] 

4A 2 = -2aLU K uj 0 /uj 2 r - (i h 0 ul/ul r )s 2 , (47) 

T + Ld^i/ujQ + 2cU'7’CJ k /(Uq — bo<-^T / ^0 T 2Xu!T^tr/^0 = 0) (48) 

and u> 2 r = u 2 r /(2/jj). Equations (l471) - (l48l) constitute the electrostatic eigenvalue equation, 
they determine A and ujq which have complex values. After writing Eq. (1171) to first order, 
we can calculate Su> so that u = oj 0 + Scu, with 8uj/ujq = 0(A)- Since the zeroth order 
operator acting on 0 (o) , £q 0) = — (r + uj^/loq) - [2u t oj k (z) - bu T /uj 0 ] + u T v 2 hi /(2l 2 ,u^)dl, 

is self-adjoint, we obtain 

dz^C^l^]) \ 


5lo 

UJq 


dzq 



U}T V thi/l\\ d 2 f 0J* e 

uo 2uj'l dz 2 V uq 


i\) 


(i) _ 


(jJT 


fkjje 

2 TTji 


(49) 




where 


r {0) _ 0 UJT V thi/ l \\ d 2 0J*i a uj t Uk{z) , u 

— o — o 77 ~ 4 o I 


-'O v* ^0 ^0 


(50) 


ujq 2ujq dz 2 

Note that the expression for 5uj only requires knowledge of the eigenfunction (/A°) to zeroth 
order. To perform the integrations in Eq. (1491) . we need the first order electromagnetic 
component, given by Eq. fl20jl . We find 


v> (1) = 


A 


b o vL/q J 0 

with 2 fi = —boS 2 ojT/{oJo\) — 2 u> k lutci/ (Acuq) 2 , and 


z /ize Xz2 + vErf fx/Az) 
dz - V J 


1 + S 2 Z 2 


(51) 


v = 


7T 


4A 


b o— (s 2 + 2A) + 2 — ( a _ 2A) 


(U 0 


UJn 


(52) 
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Thus, the electromagnetic correction to the electrostatic ITG for finite shear is 


5UJ 2 fu T , A , T , , , , A w* e , 

- — \/ — < -A— [I-IJ 3 + V Jo + T I 1-1 X 

£U 0 V VT I U 0 b 0 V ^0 




n Pi r 7 . T 1 A ^ U T , 

A A ALAn ^ 2 U 1 2 ug V 2 A 



UJT V thi/l || . 1/2 . /<^*i . ,WrW K Wt, 

3 -+ — + 4 —5-A . xl/2 

Wo 2 cu5 V Wo 025 cuq / A 1 / 2 


1 / , „2 

-- I 4——a + — 60 S 




0 


CJo 


A 3 / 2 


(53) 


with Ji = / 0 °° dzErf 2 ^\/A^ (1 + s 2 A) 1 , J 2 = / 0 °° dzzErf (VXz'j e A2|2 (l + sA 2 ) 1 , and 
Wo°° dzz 2 e 2Az2 (l + sA 2 ) 1 . We find an analytic closed form of Eq. (153|) if we introduce 
the Pade approximants for the two asymptotic limits s 2 A and s 2 -C A. For the integral 
Jij we hnd 


a 1 / 2 / 2 

2S V 7T / ' 


+ ^!og( 1 +A2) 


s 

A 1 / 2 


1 + W 

see Fig. (J7J). 

The Pade approximant of J 2 , f° r the two asymptotic limits, s 2 A and s 2 -C A, is 


(54) 


AJ 2 (s 2 ,A)~ 

see Fig. (J 8 J) . 

For the integral J 3 , we have 


1 

2P2 


1 _ 5|i 

4 A 


1 - log (1 + V2) (f) 


1/3' 


1 + — ( — ) 
1 + 8V2 V A / 


2x7/3 


(55) 


\ 3 ' 2 J 3 (s 2 ,\) ~ 


1 

3 s 2 

1- 

CO 

<M I 
<cc| 

loo I 

1 ^ 1 

32 A 

L w j 


(56) 


i + i!(f ) 7/3 ’ 

see Fig. (]9]) . 

In Figs. m~m, we pl°f the contour of SAA = 0 for a Tokamak s — a equilibrium M. 
and for a general equilibrium with arbitrary a [2lJ. We have introduced the familiar nor¬ 
malised pressure grandient parameter of ideal MHD theory, a = —R^q 2 [d'. By expanding for 
small argument the poloidal dependence of the magnetic drift frequency in an axisymmetric 
field, uj s d ymm oc cos 9 + [s# — a sin 9} sin 9, we hnd a = 0.5 — s + a. We then use /A« = 0.04, 
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v th.i/{V2l\\W*i) — 10 -2 , b 0 = 0.1, r]i = 10, r = 1, /3j = 10 -4 , and Eqs. (I54H56D . We checked 
a posteriori that for these parameters 3?[A] > 0. In Fig. (ITOl) we see that, for the s — a 
equilibium, the magnetic shear is stabilizing whereas, as expected, a has a destabilizing ef¬ 
fect. The critical a for destabilization is a growing function of the local shear. For a generic 
equilibrium [Fig. (ITT]) ], we find a critical length (for stabilization) of the extent of the unfa¬ 
vorable curvature along the field, 8d ~ /ya -1 / 2 . This quantity is a decreasing function of the 
magnetic shear. As in the s = 0 case of Sec. (13. an electron contribution is the main cause 
of stabilization. However, now the electron component of the parallel magnetic compress¬ 
ibility is subdominant, since it does not depend on the local shear, s, while J\ ~ s _1 3> 1, 
for s -C 1. The dominant term J\ is generated by the first order correction to the parallel 
component of the magnetic potential, ^ l \ calculated in Eq. (l5Tj) . which is responsible for 
electron parallel streaming, as evident from the electron solution Eq. (jffj) . 


VI. CONCLUSIONS. 

In the present work, we have revisited the problem of how curvature-driven ITG insta¬ 
bilities are affected by finite plasma pressure. As is well known, the latter affects both the 
equilibrium and the perturbed magnetic drifts of the ions, and these effects partly cancel 
each other. If the magnetic-field curvature is held constant while the electron + ion pres¬ 
sure is increased, the equilibrium VB-drift is reduced in bad-curvature regions, see Eq. (6), 
which is stabilizing. On the other hand, the hnite ion pressure gradient also introduces a 
new B x VhB|| ion drift, which is destabilizing by a mechanism identihed in Fig. 1 and 
tends to cancel the stabilizing effect of the ion pressure gradient (if the curvature k is held 
constant). There remains, however, the stabilizing action of the equilibrium electron pres¬ 
sure gradient, which stabilizes the curvature-driven ITG mode at an electron beta of order 
(3 e ~ Lt £ /R■ This scaling, heuristically derived in Sec. II, is confirmed quantitatively by 
the solution of Eq. (ITU) and shown in Figs. (EJ) and (JH]). The general dispersion rela¬ 
tion in Eq. (J23]), however, also captures the ion /?,; for the destabilization of ideal MHD 
modes, /3 MHD ~ L T ./(2q 2 R). The toroidal branch of the ITG can be completely stabilized 
f° r A A P mhd - The solution of Eq. (I23lh plotted in Fig (J2|), shows such stabilization. Fig¬ 
ures (|3j) and ( 1 j I ) confirm the scaling of the critical beta for stabilization j3^ r ^ j /!y /(2 q 2 R)- 
The comparison of /3 e ~ Lt £ /R and /3 MHD ~ Ljv/(2g 2 i?) determines which effect is more 
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important in the electromagnetic stabilization of the ITG mode. In a gyrokinetic code, this 
phenomenology is fully accounted for only if the magnetic-held perturbation 5B\\ is included. 
In particular, in its absence, the destabilizing action of the B x S75B y ion drift will be missed 
and the code will tend to underestimate curvature-driven ITG instability. 

A third critical for stabilization might be caused by the presence of a fast particle 
species. We argue that the scaling for (3J lt should be in qualitative agreement with /3^ rit ~ 
LT e /R, due to some similarities in the response of a fast population and electrons. Also 
in this case, a key role is played by the stabilizing action of the equilibrium fast particle 
pressure gradient. 

The results obtained from the local dispersion relation Eq. (I2B1) are valid when the mag¬ 
netic shear and the finite extent (along the held) of the bad-curvature region are negligible, 
unlike in a toroidal device. When these are retained, we have shown that the effect of a 
small plasma pressure gradient can be determined by perturbation theory. Since the unper¬ 
turbed (zero -(3) operator is self-adjoint, the amount of stabilization or destablization can be 
determined without calculating the perturbed eigenfunctions. The resulting expression ()53l) 
is nevertheless complicated but predicts that the extent of the unfavourable curvature along 
the magnetic held needed for electromagnetic stabilization is a decreasing function of the 
magnetic shear. 
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VII. APPENDIX 


As mentioned in the Introduction and at several places in the literature [5|, [13, [14, 23], 
the destabilizing effect of the B x S75B y drift is approximately cancelled by the stabilizing 
influence of the finite-/? modification of the equilibrium drift velocity. Mathematically, this 
cancellation can be seen directly from the kinetic equation for the distribution function 
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f = fo + Sf, where the following combination of terms appear in first order, 


v, ■ Vh/ + hv d • V/ 0 . 


Substituting the expressions (I5ll6l) for and from the Introduction gives 

2 

v d ■ VSf + <5v d ■ V/o = V K • VSf - ^|b • (Vp x VSf - V/o x VSp ± ). 

The terms within the brackets obviously have the tendency to cancel, and indeed do so 
exactly when the divergence of the current is calculated, which is effectively what is done in 
deriving Eq. (1201) . If we multiply by the charge, integrate over velocity space and sum over 
all species s, these terms disappear: 


^ e s J ( v ds ■ V5f s + Sv ds • V/. s0 ) d 3 v = ^ e s j v KS • V5f s 


d 3 v. 
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Figure 1: Physical mechanism of the curvature-driven ITG instability. In equilibrium, the ions 
drift down (in the negative y-direction) and the electrons drift up. A sinusoidal displacement £ of 
the plasma results in a positive ion pressure perturbation, 5p > 0, at k±y = (2 n + 1/2)7 r and a 
corresponding negative perturbation at k±y = (2 n — l/2)ir. Since the magnetic drift velocity is 
proportional to energy, there is an excess of ions drifting downward where 5pi > 0. Ion guiding cen¬ 
ters will therefore accumulate at k±y = 2nir and a corresponding deficit forms at k±y = (2 n + l)7r, 
which creates an upward electric field at k±y = (2 n + 1/2)7r and an E X B drift that reinforces 
the initial perturbation. An instatbility thus arises. Furthermore, pressure perturbations are an¬ 
ticorrelated with perturbations of the magnetic field strength, 6B^. The latter therefore cause a 
horizontal perturbed B X VB drift at k±_y = 2nn, which enhances the accumulation of positive 
charge at these points and strengthens the instability. 
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Figure 2: Real and imaginary part of the roots of Eq. (1231) for r = 1, fc 2 Z|| = /c 2 gi? = 0.5, b = 0.1, 
R/L Ti = 5, ut / [ythi/qR) = VbR/L Ti , co K /(v t hi/qR ) = \4>, co^r] e = cot- The unstable root at low -fy 
is the toroidal branch of the ITG mode. The mode is stabilized for /3j ~ 1.5%. 
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Figure 3: The critical /% for stabilization, as a function of ojt / (vthi, / qR) = \Jb/2qR/LT i , calculated 
from the solution of Eq. (l23l) . All parameters are as in Fig. ([2]). The line is from Eq. (1251) . 
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Figure 4: The critical /% for stabilization, as a function of k z ln = k z qR , calculated from the solution 
of Eq. (123[) . All parameters are as in Fig. (|2|). The line is from Eq. (1251) . 
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Figure 5: The critical /% for stabilization, as a function of ujt / {vthi / qR), calculated from the solution 
of Eq. (1231) . Here k z l\\ = k~qR. = 0.001, ui K /(v t hi/qR ) = 0.25, r = 1, b = 0.05, and q = 1.58. The 
line is from Eq. (12Uj) with A = 2. 
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Figure 6: The critical /% for stabilization, as a function of uj k /( vthi/qR) = Qy/b/2, calculated from 
the solution of Eq. (l23l) . Here k z qR = 0.001, UT/{v t hi/qR ) = 10, r = 1, and q = 1.58. The line is 
from Eq. (12111) with A = 2. 





Figure 7: Comparison of the numerical solution (points) and analytical solution (line) of J\ as 
defined in Eq. f)54|) . 
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Figure 8 : Comparison of the numerical solution (points) and the analytical solution (line) of J 2 as 
defined in Eq. (j55[) . 
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Figure 9: Comparison of the numerical solution (points) and the analytical solution (line) of J 3 as 
defined in Eq. (l56l) . 
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Figure 10: Contour plot of 7y[<5u;] = 0 for an s — a equilibrium (a = 0.5 — s + a). Here uj K /u*i = 
0.04, v t ,hi/(\/2l\\(jJ*i) = 10~ 2 , bo = 0.1, i)i = 10, r = 1, = 10 -4 . Toroidal electrostatic drive: 

9[wo/w«] 2> 3?[wo/w*i], and 3?[A] > 0. 
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Figure 11: Contour plot of = 0 for a general equilibrium with a = 0.5 —s+a, and a = — Roq 2 (3' 
normalised pressure grandient parameter of ideal MHD. Other parameters are as in Fig. cnD. 
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